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Abstract 

In this report, the expHcit probabiUty density functions of the random EucUdean distances associated 
with regular hexagons are given, when the two endpoints of a link are randomly distributed in the same 
hexagon, and two adjacent hexagons sharing a side, respectively. Simulation results show the accuracy 
of the obtained closed-form distance distribution functions, which are important in a wide variety of 
applied sciences and engineering fields. In particular, hexagons are often used in wireless communication 
networks such as the cellular systems. The correctness of these distance distribution functions is validated 
by a recursion and a probabilistic sum. The first two statistical moments of the random distances, and 
the polynomial fits of the density functions are also given in this report for practical uses. 

Index Terms 

Random distances; distance distribution functions; regular hexagons 

I. The Problem 

Define a "unit hexagon" as the regular hexagon with a side length of 1. Picking two points 
uniformly at random from the interior of a unit hexagon, or between two adjacent unit hexagons 
sharing a side, the goal is to obtain the explicit probability density function (PDF) of the random 
distances between these two endpoints. 

II. Distance Distributions Associated with Regular Hexagons 

A. Random Distances within a Unit Hexagon 

A unit hexagon can be decomposed into three congruent rhombuses, as shown in Fig. \T\ 
Define each one of these rhombuses as a "unit rhombus", i.e., with an acute angle of | and a 



2 




side length of 1. Here we first give the results of the random distances associated with these 
rhombuses, and then use a probabilistic sum to combine them and obtain the distribution of 
random distances for a unit hexagon. 

1) \PQ\: The probability density function of the random Euclidean distances between two 
uniformly distributed points that are both inside the same unit rhombus is [jjp| 



+ 



SlL^l rf2 _ 16 



9^3 7 

d2 



d + 



2tt 
V3 



< < 



+ f )sin--^ 



^ f 1 - C) sin 




||-(i-^)c^^ + v/4^-a-l l<ci<v/3 



-1 



otherwise 
(1) 



2) \EF\: The probability density function of the random Euclidean distances between two 
uniformly distributed points, one in each of the two adjacent unit rhombuses sharing a side but 



with different orientation, as illustrated by E and F in Fig. [H is 
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3) Final Result: If we look at the two random endpoints of a given link inside a unit hexagon 
as shown in Fig. [H they will fall into one of the two following cases: both endpoints are inside 
the same unit rhombus, with probability |; each endpoint falls into one of the two adjacent 
rhombuses sharing a side, with probability |. Therefore, given the results in Section Hl-All and 
III-A2[ the probability density function of the random Euclidean distances between these two 
endpoints is ^fo^id) + ^foiid). We hence have the following: 

[Random distances within a unit hexagon] The probability density function of the random 
Euclidean distances between two uniformly distributed points that are both inside the same unit 
hexagon is 
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Fig. 2: Random Points between Two Adjacent Hexagons: Different Cases with Rhombuses. 



The corresponding cumulative distribution function (CDF) is 
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5. Random Distances between Two Adjacent Unit Hexagons Sharing a Side 

Given two unit hexagons that are adjacent to each other, as shown in Fig. [2l they can be 
decomposed into six congruent rhombuses in a similar way as that in Fig. [T] Here we simply 
give the results for the random distances between these unit rhombuses with difference placement, 
as shown in Fig. |2l At the end we combine them by a probabilistic sum, in order to obtain the 
distribution for the random distances between two adjacent hexagons sharing a side. 
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1) \DE\ 
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5) \HE\: 
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6) Final Result: Given the resuhs in Section Ill-All to III-B51 when the two endpoints of 
a given link fall into one of the two adjacent hexagons sharing a side, the probability density 
function of the random distances between these two endpoints is | [/Dal*^) + fo^d) + /^^((i)] + 
§ [foiid) + foiid) + foaid)], using the similar reasoning as that in Section III- A3 1 With this 
probabilistic sum, we thus have the following: 

[Random distances between two unit hexagons] The probability density function of the 
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random Euclidean distances between two uniformly distributed points, one in each 
adjacent unit hexagons sharing a side, is 
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The corresponding CDF is 
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Note that although unit rhombuses and unit hexagons are assumed throughout (fTI)-(fTTI). the 
distance distribution functions can be easily scaled by a nonzero scalar, for rhombuses or 
hexagons of arbitrary side length. For example, let the side length of a regular hexagon be 
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Fig. 3: Distributions of Random Distances Associated with Hexagons. 



s > 0, then 

FsD{d) = P{sD <d)= P{D <-) = Fn{-). 

s s 

Therefore, 

Uid) = F'^i-) = (12) 

s s s 

where foi-) can be © or (flOl) . 

III. Verification and Validation 

A. Verification by Simulation 

Figure [3] plots the probability density functions of the two random distance cases given in 
([3]) and (fTOl) . respectively. Figure |4] shows the comparison between the cumulative distribution 
functions (CDFs) of the random distances, and the simulation results by generating 2, 000 pairs 
of random points with the corresponding geometric locations. Figure |4] demonstrates that our 
distance distribution functions are very accurate when compared with the simulation results. 



12 



1 1 l^>H-rii' 1 


V 


/ 


J 

/ 

/ 






/ '''' 

i ^ 
r ^ 




4 ^ 

f 

J «' 
/ «' 

- / ^ 


Within Hex Distribution 

+ Within Hex Simulation 

Between Hex Distribution 

" Between Hex Simulation 







0.5 1 1.5 2 2.5 3 3.5 4 

Distance - d 



Fig. 4: Distribution and Simulation Results For Random Distances Associated with Hexagons. 




Fig. 5: Partial Recursion through Hexagons and Rhombuses. 



B. Validation by Recursion 

As shown in Fig. [5l a hexagon with a side length of 2 can be decomposed into three small 
hexagons with a side length of 1, and three rhombuses A, J and K, each with a side length 
of 1 as well. With the scale transform in (fT2l) . the distance distribution in the large hexagon is 
f2D{d) = |/DHj(i)- On the other hand, if we look at the two random endpoints of a given link 
inside the large hexagon, they will fall into one of the three following cases: i) both endpoints 
fall inside one of the three small hexagons (BDF, ECG or HIL), with probability | x |; ii) one 
of the endpoints falls into one of the small hexagons, and the other endpoint falls into one of 
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the three rhombuses A, J or K, with probability 2 x | x i; iii) both endpoints fall into one of 
the rhombuses, with probability j x j- 

Each of these three cases includes several more detailed sub-cases as follows: 
Case i) Given the location of the first endpoint of a particular link, the second endpoint will fall 
in the same hexagon as the first one with probability |, and in one of the two adjacent 
hexagons with probability |. The unconditional probability of these two sub-cases are 
^x| = ^and^x| = |, respectively. 
Case ii) Without loss of generality, suppose the first endpoint is located in A. By symmetry, the 
second endpoint falls into any one of the four rhombuses B, D, F and H with probability 
|, and into L with probability |. Thus the unconditional probabilities for \AB\, \AD\, 
\AF\ or \AH\ are all | x | = ^, and | x | = ^ for \AL\. 
Case iii) If the first endpoint is in A, then by symmetry, the second endpoint is still located in 
A with probability |, and in either one of J or K with probability |. The unconditional 
probability of these two sub-cases are ^ x | = ^ and x | = ^» respectively. 
In short, we have the probabilistic sum as 

f2D{d) = ^fD^M) + lfl^nAd) + ^[fnM+fDM+fDM+fDe{d)] 

+ Y^[fDM) + fDM)] + ^fDM), (13) 

where fosid) is the density of \AL\ in Fig. [51 which is the only distance distribution function 
that has not been given yet. We will do so in the immediate following. 
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1) \AL\: 
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2) Validation: In order to confirm that the two definitions of f2D{d) at the beginning of Sec- 
tion |in:B] are equivalent, i.e., ^/dhj(c?) + Ifoujd) + [JdM + foM + foM + JdM] + 
^ [f]:,^{d) + fDsid)] + ^foiid) is equal to ^fon^i^), we verify them mathematically as follows. 
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while the probability density functions are for all other cases. Thus, 

f2D{d) = ^fD,,{d) + lfn,^{d) + ^[foM + fDM] + ^fDM 
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+ 1 -— - - - 4Vd2 - 3 + — V4d2 - 27 

.9^3 3y 6 3^3 2 



V3V3 V3/ d 



d^ 1 \ . ,V3 5^3 . _i 3v^ 

— ;= H ^ sm — - H sm — — 

3^3 2V3J 2d 2 2d 



■-^^ — Vd^-3+ -V4d2 - 27 p - - 

4 6 2v^ 4 
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JdM) 



d 
12 

+ 



1 f2d^ . 1 V3 4 X /g 

^ — + 1 sin"^ ^ — + 2 sin"^ -V 

^V3 J 2d J d 



(3V3 



2^2 /-\ . _i 3v^ , v/4rf^ 
+ 3v3 I sm ' 



2d 



4:Vd^^ + ^ V4d2 _ 27 - ^ 
2 3^3 



12 



2 /(i^ o\ . 1^3 4 /ci2 X . . 3v^ /I 27r \ 
^ — + 8 sin"^ ^ — + 12 sin"^ + - + — 7= d^ 



+6 - -VM^ -27 + ^ + 7 
3 3^3 
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Thus, 



f2D{d) = lfD^,id) + ^fnM) + ^[fDM + fnsid)] 



d 
6 



27T fd 



2{d/2) 3^3 V2 



ix) 2\/3 < d < \/l3 : 



12 



2 

V3 V3 



I — + 12 sin"^ -f- + sin-i 



19 



16 



2d d 
16n 25 



2 47r 
3^ 9^3 



12 



fDeid) 



d 
6 



(3V/3 ' V3 



sm 



2x/3 _ / 



d 



d^ 1 \ . _i \/3 5\/3 . _i 3\/3 

\ — ^ ^ ^ sm — — H — sm — — - 

V3v^ 273/ 2d 2 2d 



1 TT \ ,2 V4d2 - 3 5 r— — — „ r-r^ — — 3l7r 9 

- + ^)d^-- — + -V4d2-27 + 2Vd2-12 p-- 

6 9\/3/ 4 6 6^3 4 



d 
12 



1 /2d2 \ . ^ 



2d2 



1 27r 
3^ 9V3 



2d V3V3 
\/4d2 - 3 3 



. . -i3v^^ 8 . _i2v^ 
3V3 sm —— + — sm — — 
2d d 



+ ^V4d^+^Vd^-i^ 



Y^^sid) 



_d 
12 



4 ^t^' 



.1 3\/3 _ 2 fd^ 
'2d~^\J 



+ 8 sin" 



-1 2x/3^/ ^^^^2 



d 



3^3 



-^V4d2 - 27 - 4v/d2^+ + 17 
3 3^3 
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Thus, 

f2Did) 



d 
6 

d 
6 



1 



TT 



10 



IGtt 



-4 



/ 47r 



20 

y 



3- 



IGtt 
3^3 



x) V13 < < 4 

1 
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Thus, 

f2Did) 



d 

12 



h 16 sm — — 

V3 v3 ; d 



1 2tt 
3 ^97! 



3 



12 



327r 
3^3 



7! 



d 

6 



Til 3 



- 16j sin ^ 

.1 73 _ 
d/2 



2V3 / 



TT 



47r 2 

V!^ 3 



V9v^ 6 
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2 



3 



12 



IGtt 
3^3 

IGtt 
3v^ 



In summary, we have f2D{d) = |/DHi(i) by recursion, and the probabilistic sum ^/dhjI*^) + 
lfD^^{d) + ^ [fD,{d) + + /d3(^^) + fnAd)] + M t/^rl^^) + foAd)] + ^/dhj(^^) is equal 

to |/dhj (f ) in all the cases discussed above. The results are a strong validation of the correctness 
of the distance distribution functions that we have derived. 



IV. Practical Results 



A. Statistical Moments of Random Distances 

The distance distribution functions given in Section |II] can conveniently lead to all the statistical 
moments of the random distances associated with hexagons. Given fDu^id) in ([3]), for example. 
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TABLE I: Moments and Variance — Numerical vs Simulation Results 



Geometry 


PDF/Sim 






VarD 


Within a 
Single Hexagon 




0.8262542775s 


0.8333333333s 


0.1506291100s^ 


Sim 


0.8263306317s 


0.8335924725s 


0.1507701596s^ 


Between two 
Adjacent Hexagons 




1.8564318344s 


3.832947195s 


0.3866080394s^ 


Sim 


1.8583366966s 


3.8326819696s 


0.3792666917s^ 



the first moment (mean) of d, or the average distance within a regular hexagon, is 



M 



(1) 



^fou^ {x)dx 



7v^_ 7 1 
~W ^ 90 ^ 60 



28 In ( 2^3 + 3 ) + 29 In ( 2^3 - 3 



and the second raw moment is 



M 



(2) 



from which the variance, or the second central moment, can be derived as 



M 



(1) 



1 2 



0.1506291100. 



0.8262542775, 



When the side length of a hexagon is scaled by s, the corresponding first two statistical 
moments given above then become 

5s 



M 



(1) 



0.8262542775s, M. 



(2) 



D 



Hi 6 



and VarDt,^ = 0.1506291100s^ (15) 



Table U lists the first two moments, and the variance of the random distances in the two cases 
given in Section UIl and the corresponding simulation results for verification purposes. 



B. Polynomial Fits of Random Distances 

Table HI] lists the coefficients of the high-order polynomial fits of the original PDFs given in 
Section [III from the highest degree (degree-10 for ([3]) and degree-20 for (flOl) ) to rf", together 
with the corresponding norm of residuals. Figure [6] (a)-(b) plot the polynomials listed in Table HI] 
with the original PDFs. From the figure, it can be seen that both polynomials match closely with 



25 



TABLE n: Coefficients of the Polynomial Fit and the Norm of Residuals (NR) 



PDF 


Degree 


Polynomial Coefficients 


NR 




10 


10^ X [-0.0146710 0.136604 - 0.538052 1.167903 
—1.525478 1.230615 —0.605940 0.175147 — 0.043772 

0.025830 - 0.000025] 


0.075608 




20 


10^ X [0.00000035 - 0.000013 0.000207 - 0.002094 
0.014469 - 0.072522 0.272508 - 0.782682 1.736254 
-2.986406 3.976655 -4.072372 3.169347 - 1.841066 
0.778001 - 0.230634 0.045522 - 0.005534 0.000394 
-0.0000103 0.00000007092] 


0.191157 




the original PDFs. These high-order polynomials facilitate further manipulations of the distance 
distribution functions, with a high accuracy. 

V. Conclusions 

In this report, we gave the closed-form probability density functions of the random distances 
associated with hexagons. The correctness of the obtained results has been verified by a recursion 
and a probabilistic sum, in addition to simulation. The first two statistical moments, and the 
polynomial fits of the density functions are also given for practical uses. 
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